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ABSTRACT 

We suggest that the Baldwin effect is a result of the spectral dependence of the 
line-driving ionizing continuum on the black hole mass. We derive a relationship 
between the mass of the central black hole and the broad emission-line luminosity 
in active galactic nuclei (AGN). Assuming the UV spectrum of AGN is emitted 
from an optically thick medium we find an expression for the characteristic energy 
of the "UV bump" in terms of the observable luminosity and emission-line width. 
We show empirically and analytically that the bump energy is anti-correlated with 
the black-hole mass and with the continuum luminosity. Our model reproduces the 
observed inverse correlation between equivalent width and continuum luminosity, 
yielding an explanation of the Baldwin effect from first principles. The model gives 
a good fit to the Baldwin effect of the CIV line for a mean quasar EUV spectrum 
(Zheng et aZ.1997) and for several model spectra. The model also predicts a correlation 
between the strength of the Baldwin effect (the slope of the equivalent width as a 
function of luminosity) and the ionization potential, consistent with recent data. 



Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — quasars: 
general — black hole physics — accretion disks 



1. INTRODUCTION 

Baldwin (1977) showed that the equivalent width of the CIV emission- line decreases with 
increasing continuum luminosity, a correlation known as the 'Baldwin effect'. This result has 
been confirmed for a large set of AGN and a wide luminosity range (Kinney, Rivolo and Koratkar 
1990). A similar relation has been found also for other broad emission- line s. Several models have 
been proposed to explain the Baldwin effect , such as a decrease of the ionization parameter or of 
the covering factor with luminosity (Wu, Boggess and Gull 1983) or the accretion disk inclination 
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(Netzer 1985). Nevertheless, the origin of the Baldwin effect is still not understood (for a recent 
review see Shields and Osmer, 1998). 

The widest accepted explanation for the energy source of quasars and Active Galactic Nuclei 
(AGN) is the accretion of matter onto a supermassive black hole in the center of the host galaxy. 
The central mass implied by a variety of observations (see Wandel 1998) is of the order of 
10 6 — 1O 9 M0. The angular momentum directs the accreted material to form a disk around the 
central black hole and the properties of such an accretion disk may be calculated under a few 
basic assumptions (e.g. Sakura and Sunyaev, 1973). In particular one can calculate the luminosity 
and spectral distribution of the radiation from the accretion disk, given the black hole mass and 
the accretion rate. 

The thermal emission from the accretion disk around a central mass of the size implied 
for AGN peaks in the UV, and ionizes the broad emission-line gas. The precise shape of the 
continuum spectrum from an accretion disk depends (among other factors) on the central mass. 
Since the continuum radiation ionizes the line-emitting material and drives the broad AGN 
emission-lines , the emission-line strengths depend on the continuum spectrum and therefore on 
the central mass. 

In particular the spectral shape of the ionizing radiation can determine the equivalent 
width (the ratio between the energy flux in the line to the continuum energy flux at the line 
frequency) of the emission- line s. 

We suggest that the Baldwin effect is a result of the spectral dependence of the line-driving 
ionizing continuum on the black hole mass. In the next section we present constraints on 
the black hole mass; section 3 describes the temperature and spectrum of the basic thin 
accretion disk model. Section 4 gives an empirical relation between the model parameters (M 
and M) and the observed spectrum and luminosity, while section 5 relates these parameters to 
the calculated spectrum for various accretion disk models. In section 6 we derive the dependence 
of the line luminosity on the physical conditions in the gas. In section in 7 we compose all these 
elements to produce a relation between the equivalent width and the mass, and in section 8 we 
show that the model predictions are indeed consistent with the observed Baldwin effect for CIV 
and other lines. 

2. Basic relations for accreting Black Holes 

Any compact, accretion powered radiation source must obey several fundamental relations: 

a. The Eddington limit: in order to maintain steady spherical accretion the luminosity 
must be less than the Eddington luminosity, 

L < L Edd = ATiGMm p c/a T = 1.310 46 M 8 erg/sec 
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or 

M 8 = 0.7^-^46 (1) 
where M 8 = M/1O 8 M , i] = L/LEdd is the Eddington ratio, and L 46 = L/10 46 erg/sec. 

b. The black body temperature: if a luminosity L comes from a region of radius R and 
temperature T, then 

L < 4ttR 2 <tT*. (2) 

If a spectral feature at a photon energy of E is due to black body emission from a region of 
size R, then the temperature is given by E ~ 3kT, and E ~ (10 eV)(T/W 5 K), so that 

£ ~< (10 eV)L\( i A M 8 1/2 (R/5R s )-^ 2 . (3) 

and reversing the relation above we have an lower limit on the black hole mass (for a true 
black body spectrum this becomes an approximate equality): 

M 8 ~> 2E R 2 d L 4e 1 / 2 (R/5R s )- 1 (4) 

where R s = 2GM/c 2 « 3 10 13 M 8 cm is the Schwartzshild radius and E^yd is the spectral 
energy in Rydbergs. (Note that L^q is the bolometric luminosity; if the observed (say optical) 
luminosity is used, a bolometric correction must be included). 

Combining eqs. [I] and [| we can eliminate M obtaining (for black body emission) an estimate 
of the size: 

§- « 13L- b 1/2 V E Ryd - 2 (5) 



R 



s 



c. Variability: the shortest time scale for global variations in the luminosity is the light 
travel time across the Schwartzshild radius, hence if the luminosity is observed to vary significantly 
on a time scale St, the black hole mass has to be 

Ms < TrJf~ > (6) 
lO-Vsec 

where r is the effective radius of emission in units of R s . 

However, for continuum changes near the feature of spectral energy E the light travel time is 
(from eq. IB]) 

R 

Sti& — k, 1.5L^ 2 £;^" 2 days, (7) 



and the dynamical time is 



St d « ^ « 20L^y/ 2 £ i? j /( T 3 days, (8) 
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which gives time scales of the order of the observed UV variability in Seyferts and quasars. 



3. The thin accretion-disk model 

Many authors have tried to fit accretion disk spectra to the observed AGN continuum, 
implying the accretion disk parameters (M, M, a) and the black hole spin (e.g. Schwartzshild or 
Kerr) which give the best fit to the observed optical-UV continuum (e.g. Wandel and Petrosian 
1988; Sun and Malkan 1989; Laor 1990). 



3.1. The accretion disk spectrum 

The radiative energy output of the accretion disk is dominated by the release of gravitational 
potential energy, the rate of which is GMM/r. A more accurate calculation yields for the flux 
emitted per unit area 



3GMM 
L{R) ~ ~8ttR^~ 



R X 1 / 2 
~R 



(9) 



where Ri n is the inner disk radius. In the outer part of the thin disk the opacity is dominated 
by true absorption and the local spectrum is a black body spectrum. For this part of the disk, 
comparing eq. || to the black body radiation flux gives for the disk surface temperature 

T w{-^w) ~ 610 U) r 1 K < 10 > 

where r = R/R s and 

m = Mj Mem « 2(M/M Q yr~ 1 )M 8 1 (e/0.iy 1 

is the accretion rate in units of the Eddington accretion rate , rh = 1 at the Eddington limit. 
When the accretion rate approaches the Eddington rate, the intermediate disk region becomes 
dominated by electron-scattering, the spectral function will be of a modified black body and the 



surface temperature is higher than given by eq. 1C. At still smaller radii, the pressure in the disk is 



dominated by radiation pressure, rather than by gas pressure, and the thin disk solution becomes 
thermally unstable. In that inner region the thin disk solution is probably not valid, and has to be 
replaced by a hot disk solution (e.g. Wandel and Liang 1991). It turns out that the intermediate 
modified black body region is relatively narrow, and close to the black body - electron scattering 
boundary the spectrum is nearly a black body one, so both regimes may be approximated by the 
black body solution. 

The spectrum of a multi- black body accretion disk is given by integrating over the entire 

disk, 

L v rs f ° Ut 2irRB u \T(R)]dR (11) 
JRt 
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where B V {T) is the Planck function and Rt is the transition radius from the intermediate to 
the inner radiation-pressure dominated region (for m < 0.02 the black body region extends down 
to the inner edge of the disk and R t = Ri n ) . 

Since the B V (T) has a sharp peak at hv co ~ 3kT and cuts off at higher frequencies, the highest 
frequency of the black body part of the disk spectrum comes from the radius Rt, with the highest 



temperature for which the disk is still optically thick. Eq. 11 gives a spectrum which depends on 
the radial extent of the black body part of the disk. If that part is extended (R ou t/Rt » 1) the 
spectrum is almost flat (~ v 1 / 3 ) and cuts off beyond hv co « 3kT(R t ) (or 3kT{5R s )/h for m < 0.1). 
If Rout I ' Rt ~ a few, the spectrum will be merely a somewhat broadened Planck spectrum, but if 
the accretion disk is extended over a large radial range there is a significant flat part. In that case 
a better approximation than the Planck spectrum is 

I/\V3 / y 



Lv ~ A — exp[ (12) 

where A and v co are the normalization and cutoff frequency For a Kerr black hole Malkan (1990) 
finds v w = (2.9 x 10 15 i?z)M 1/ 1 4 M 8 ~ 1/2 , where M .i = M/0.1M Q yr- 1 . This can be written as 



hv co = (6 eTOm^M- 174 = (20 e^L^M" 172 (13) 

(note that here L^q is the observed luminosity, and we have included a bolometric correction of 10. 
Also note that for a Kerr black hole the efficiency in the definition of rh is e = 0.4). The relation 
between the cutoff frequency and the inner disk temperature in the two cases of face-on Kerr and 
the black body approximation is 

hv co = OMkTmax = 2.5kTsB 

(Malkan 1990). 

Although the soft X-ray excess suggest a higher thermal temperature in some objects, In 
general the characteristic temperature is likely to be related to the UV spectrum; While in 
the UV spectrum is produced by thermal emission from an optically thick gas (with eventual 
reproccessing), the X-rays are more likely to be produced by a hotter medium due to processes 
such as Comptonization, a two-temperature disk (Wandel and Liang 1991) or hot corona (Haardt, 
Maraschi and Ghisellilni 1994; Czerny, Witt and Zycki 1996). 



4. Deriving M and M from the UV spectrum 

We want to find an approximate relation between the black hole mass and accretion rate and 
the continuum spectrum. The accretion rate is determined straightforwardly by the luminosity, via 
the relation Lj on = efi on Mc 2 , where fi on is the fraction of the ionizing out of the total luminosity. 
To estimate the black hole mass we may use the multi-black body accretion disk -spectrum. As 
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discussed in the previous subsection the UV bump cutoff frequency is determined by the highest 
surface temperature in the black body part of the disk. If the black body region extends down to 
a radius R t , and the EUV cutoff energy is E co , then eq. |] gives (we assume the inner region is 
optically thin and much hotter, so its contribution in the frequencies of interest can be neglected) 

M 8 « 2.5(E co /WeV)- 2 Ll / 6 2 (R t /5R s )- 1 . (14) 

In several accretion disk models the black body regime extends close to the inner disk edge. 
This is the case in the "/3" disk model with viscosity proportional to the gas pressure (y ~ PP ga s)- 
For low accretion rates (m < 0.1) this is true also for the "a" model with viscosity proportional to 
the total pressure [y ~ a(P gas + P ra d))- in this case, since the disk emissivity peaks at R sa 5R$ 
(for a non-rotating black hole ) we may use eq. || in order to find the cutoff energy. 

E ~> (17 eV)L% 4 M s 1/2 Schwartzshild (15) 

(where we have introduced a bolometric correction L op t = O.IL^). 

The equivalent result for a Kerr black hole (eq. |l3|) is very similar: 

E ~> (20 eV)L l J b 4 M~ 1/2 Kerr. 

Plotting the cutoff frequency versus the luminosity one can infer the mass, and vice versa, if 
the mass can be estimated independently, it is possible to estimate the cutoff frequency (fig. 1). 
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Figure 1. The cutoff energy (for an accretion disk spectrum vs. the monochromatic luminosity 
at 5100A. Diagonal dashed lines indicate constant black hole masses, and diagonal dotted lines 
indicate the Eddington ratio. The data shown consists of PG quasars (diamonds; from Boroson 
and Green, 1994), Seyfert 1 galaxies (triangles and squares, the latter indicate AGN with BLR 
reverberation data from Wandel, Peterson and Malkan 1999, with luminosities modified for 
H = 50) and Narrow Line Seyfert Is (Boiler, Brandt and Fink, 1996) 

It is possible to estimate the black hole mass from the broad emission line profile, assuming 
the velocity width is induced by a Keplerian velocity dispersion. When the broad line region size 
is estimated by reverberation mapping, this technique is particularly reliable. Wandel, Peterson 
and Malkan (1999) have used a sample of AGN with reverberation data to calibrate the mass 
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estimate obtained from photoionization models, finding the relation 



M 8 « 0.4 ( ) 1/2 v% (16) 



\Un 



10 



where t7 is the ionization parameter, nio is the density in units of 10 10 cm -3 and vs is the H/3 
FWHM in 10 3 km s _1 . Using this calibrated relation makes it possible to estimate black hole mass 
estimates for large samples, even without reverberation data (Wandel, Malkan and Peterson, in 
preparation). Combining equations [lf| and 13 or 15 we have 



E co « (40 eV)(Un 10 )- 1/4 Vs 1 (17) 
Note that the luminosity dependence cancels out, and the dependence on the unknown parameters 



U and n is weak. Using eq. 17 and the average fiducial value of Unxo=l we find the mass and 
cutoff energy distribution for a large sample of quasars, Seyfert 1 galaxies and NLSls (fig 1). Note 
that E co is anticorrelated with the mass, and in particular that different AGN categories group 
in different regions in the L — E co plane: quasars have more massive black holes and low cutoff 
energies, Seyfert 1 galaxies have intermediate black hole masses (10 7 - 10 8 M Q ), and NLSls have 
low black hole masses (10 6 - 1O 7 M ) and high cutoff energies. This is consistent with the large 
soft X-ray excesses observed for many of the NLSls (Boiler, Brandt and Fink, 1996). 

The diagonal dotted lines in fig. 1 give the L/LEdd ratio - quasars are in the 0.01-1 range, 
Seyferts in the 0.001-0.1, and NLSls are all near Eddington , L/Lem ~ 0.1 — 1. 

A similar analyses for the Schwartzshild accretion disk has been performed by Wandel and 
Petrosian (1988), who calculated for each pair of accretion disk parameters M and M the actual 
observables, e.g. the UV luminosity and spectral index. 



5. Relating M and L to the Ionizing Spectrum 

5.1. Temperature dependence of the Ionizing Spectrum 

The ionizing power Q is essentially the photon flux above the ionization potential Ei on of the 
line in question - e.g. lRyd (13.6 eV) for Lya and 48eV for CIV. 

We shall now relate the ionizing flux (defined by Q cx Li on / < E >) to the temperature of 
the continuum-emitting gas, assuming the emission is thermal (which is certainly true in the thin 
accretion disk model, and supported by the presence of the blue bump feature in the spectrum of 
most AGNs. If E co < Ei on this frequency is on the Wien part, and the flux is sensitively influenced 
by the effective disk temperature, T e ff. In the parameter range of interest we may therefore 
assume that the equivalent width is strongly correlated with T e ff. Quantitatively, the correlation 
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is related to the logarithmic derivative 

d In B v 



(18) 



d In T e x — 1 

where x = hv/kT. For kT <1 Ryd, which is true in a wide range of cases, as we show below, eq. 



18 



gives 

d *" B » « -hv/kT, 
d In T 

namely a strong correlation with T. 

In order to relate the UV ionizing spectrum near lRyd to the accretion parameters, it is 
necessary to determine the effective black body temperature of the emitting region, which we do 
in the next section. 

Once this temperature is determined, we can estimate the ratio of fluxes in any two 
desired energy bands, in particular the flux in the ionizing energy band to the energy Ej of the 
emission-line in question, which is related the equivalent width , since the line luminosity is related 
to the ionizing flux. From eq. |l^ we have 

L vl0 n/L uI « {E I /E ion ) l ^exp[{E I - E ion )/E ]. (19) 



5.2. The Variability-Effective Temperature Relation 

We have already derived a quite general temperature-radius relation for black body emission 
from an accretion flow (eq. We can apply this relation to determine the temperature or 
radiation peak energy as a function of the observed luminosity and estimated mass. 

In order to determine a characteristic temperature we need either indicate the size of 
the emitting region (and use eq. |2|) or assume an emission model, as done below for the 
accretion disk model. The size of the emitting region can be bounded by variability analyses. 
From eqs. [2| and ^ we have 

L < 47rac 2 (5t) 2 T 4 



which gives 



E eff > WL^iUayy^eV (20) 



where td ay is the variability scale in days. X-ray data imply a linear relation between luminosity 
and variability time in AGN (Barr and Mushotzky, 1986, Ulrich, Maraschi and Urry 1997). 
Recently a time scale of the order of 30 days has been found in the X-ray PDS of NGC 3516 
(Edelson and Nandra 1999), but the fastest variability observed in AGN is of the order of hours. 



Assuming a relation, t& av ~ tlL^q, eq. 20 gives 



E eff > (10 eF)r L 1/2 L 46 - 1/4 . (21) 
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For a linear M — L relation L^q ss (ij/0.1)Mg (e.g. Wandel and Yahil 1985; Wandel, Peterson and 
Malkan 1999; fig. 1 above) eq. |2l| gives 

E eff « nr£ 1 %/0.1)- 1 ^M 8 - 1 ^eV, (22) 

similar to the result obtained from accretion disk models below. 



5.3. Effective Temperature Models for an accretion disk 

In the framework of the thin accretion disk model, we consider four models for determining 
the effective temperature. 

a. The inner disk black body temperature 

When the disk is nearly black body up to the inner edge, we can approximate the spectrum 
by black body emission near the maximum-emissivity radius, which gives 

Emax « (3 eV)kT BB (5R s ) = W(m/M 8 ) 1/4 . (23) 

As discussed in sec. 3 above, this model applies to the disk and to low accretion rates 
{fa < 0.1) in the a disk. A similar expression (eq. [i~3| ) applies to the Kerr disk. 

b. The inner disk radiation-pressure temperature 

This model applies to high accretion rates (rh > 0.1) in the a disk, provided the inner disk 
maintains the temperature profile given by the thin-disk solution. The temperature in the inner, 
radiation dominated region is given by (Sakura and Sunyaev 1973) 

T rad (r) «(2x 10 5 /f)(«^8) _1/4 r- 3 / 8 [l - (6/r) 1 / 2 ] 1 / 4 . (24) 



Substituting the maximal-emissivity radius r = 5 gives 



Emax « (4 ey)(aM 8 )- 1 / 4 (25) 

and a higher value for a Kerr black hole . The gas-radiation pressure boundary - taking as T max 
the temperature at the inner boundary between the black body and the radiation-dominated disk 

c. The black body -radiation boundary temperature 

In this case the effective temperature is taken at the boundary between the intermediate and 
inner disk regions, r m j ~ 10(aMg)°' 1 m ' 76 . 



The energy of the peak in the spectrum is then derived by substituting r m j into eq. |24j, which 
gives 

E max 2(amM 8 )- a3 eV. (26) 
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This model is adequate in the case of a high accretion rate - a disk when the inner disk does 
not follow the thin disk model but establishes a stable hot solution (e.g. Wandel and Liang 1991). 

d. The "photospheric" temperature 

In this model we consider as the disk effective temperature the temperature at the radius 
where the disk becomes effectively optically thin. The effective optical depth in the disk is given 
by 

T* = hp(K es K a ) / , 

where h and p are the disk vertical scale-height and density, respectively, and K es and K a are the 
electron-scattering and absorption opacities. Equating t* to unity and eliminating T and r by 
using the temperature profile T(r) gives 

E max « 2a- a4 (mM 8 )-°- 26 eV. (27) 

Note that the last two models give very similar results. 



6. The Line Emissivity 

The equivalent width of an emission-line is actually the ratio between the luminosity in the 
line and the continuum luminosity at the emission-line frequency. 

The line luminosity is determined by the photoionizing flux, by the covering factor of the 
line-emitting gas and by the line emissivity which is a function of the flux of ionizing photons and 
the density. The differential line luminosity emitted from a thin shell can be written as 

df 

d Lune = L ion —j Une (r)dr (28) 
ar 

where df/dr is the differential covering factor, and jiine(Q,n), the line relative surface emissivity 
(or reproccessing efficiency), is a function of the flux of ionizing photons (Q = Lj on /47rr 2 c < E >), 
the density (n). < E > is the weighted mean energy per ionizing photon. The line luminosity is 
then found by integrating eq. ^ over radius 

Liine = L ion J -^JUne [Q(r),n(r)] dr (29) 

In order to relate the line luminosity to the ionizing spectrum, we need first to express the 
line luminosity in terms of the ionizing luminosity and spectrum. The emissivity of several lines 
has been charted for large portion of the Q — n (flux and density) plane using the CLOUDY 
photoionization code in a slab geometry (Baldwin et a/.1994, Korista et a/.1995). These authors 
find that some lines, in particular CIVA1549 and OVI A1094 have a particularly simple emissivity 
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structure, with the emissivity peaking strongly (and being approximately constant) along the 
diagonal Q/n =const., that is, for the ionization parameter (defined by U = Q/n) being constant, 
U » U Q . 

Approximating 

j(Q,n)^j 5(U -Uo) 
(the locally maximally emitting model) the line luminosity can then be written as 



Lline ~ Li m f e ffjo, (30) 

f e ff is the effective covering factor, given approximately by integrating df/dr convolved with j(r) 
over the emitting layer, and the effective radius (e.g. for calculating f e ff is given by the condition 
U(r eff ) = U ,ot 

^x Hre// )]V^( _^_ )' /2 . (31 ) 



7. Relating the equivalent width to M and L 

From the definition of the equivalent width we have EWj « Lj/L c , where I denotes the line 
in question and L c is the continuum flux at the frequency corresponding to that line. In section 
5 we have seen that the thermal UV spectrum from the accretion disk has a strong dependence 
on the effective surface temperature of the disk, and that for a variety of models the functional 
dependence of the effective photon energy on the black hole mass is approximately 

E eff (M) « (3 - 6 eV)M s 1/4 (32) 

where the range in the coefficient accounts for the uncertainty in the parameters m and a (for the 
accretion disk models) or r\ and tl (for the variability-temperature relation). 



7.1. The Temperature-Luminosity Relation 

A more direct relation can be obtained with the luminosity (which is what we actually need 
in the context of the Baldwin effect ). 

The variability-temperature relation gives E e ff ~ (10 eV){TL / lday)~ l l 2 L^ 1 ^ 4 , For the 
accretion disk models we express M and rh in terms of L. Since, from the definition of m we have 
mM 8 = 1.3 x 10 46 L 6oii46 erg/sec ps (10 - 20)L optA6 , so that 

(mMg) 1 / 4 « 2L% A 
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The inner black body model (eq. |2^) gives 

E eff « (4 eVO(m/0.1) 1 / 2 ( a /0.1)- 1 / 4 L« 1/4 
the radiation-pressure model (eq. |25|) - 

E eff « (8 e^Cm/O.^^Ca/O.l)- 1 / 4 ^ 1 / 4 
two last accretion disk models (eqs. [2(| 27) can be combined into a single expression 



£ e// « (2 eF)a- 1 / 3 m- 1 / 4 M 8 ~ 1/4 

which gives 

£ e// « (3 eF)(a/0.1)- 1 / 3 L 4 - 1/4 . 

Note that the accretion disk models are not including Comptonization by a hot gas, which is 
likely to increase E e ff, so they are probably underestimating E e ff. 

Combining the variability and accretion disk estimates we can write 

E eff (L) * E C L~ 6 1/4 (33) 

with E c ~ (10 — 20)eV, allowing for a range in the parameters a, tl and m-. This theoretical result 
is consistent with the empirical distribution of E co independently found from the BLR data (fig. 
1), and with direct observations of the far UV spectrum of AGN, which show that the SEDs of 
several AGN observed with IUE and HUT peaks at 1000A(Zheng et aZ.1997, Kriss et al. 1999). 



7.2. Power Law Continuum 

The observed SEDs suggest that the decrease in the far UV is a power law (~ i/( _ ( 15-2 ) 
rather than exponential. Assuming a UV spectrum of the form 

L v ~ u^ 1 for hu < Eb 

L u ~ v~ a for hu > Ef) 

and Ei < E b < E ion gives 

EWi cx L ion /vL u (Ej) « {E b /E ion ) a ~ l 
because the line luminosity is proportional to the ionizing luminosity Li on ~ vL v {Ei on ). 
Taking E b = E e ft we then have 

EWi cx L im /Li « (EiojEt)- 1 « {E c lE ion )- x L~ {a -^l\ (34) 
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For a ~ 1.5 - 2 (Zheng el a/.1997) this gives EW oc L~ 7 with 7 ~ 0.13 - 0.25. 

While for a pure power law, this would imply the same slope for all lines, if the spectrum 
is curving to steeper slopes at higher photon energies (as could be the case for a partially 
Comptonized Wien continuum), lines with higher ionization potential (e.g. OVI) would have a 
steeper slope and vice versa - also consitent with the data. 



7.3. Wien Continuum 



From eq. 12 



EWi = L uion /L uI j f e ff oc {E I /E ion ) 1 ^exp[{E I - E ion )/E eff ) (35) 
Combining these two equations we have 

EWj oc exp ( El (36) 

Subject to these uncertainties, we can try to estimate the Baldwin effect quantitatively. 
For the CIV A1546 line (£/ - E ion ) = (8 - 48) eV = -40 eV, which gives 

../aasn/EWciv = EW0 C ivexp(-40 eVE~ l L% A ). (37) 

The coefficient can be determined by normalizing the model to the low luminosity en to the data, 
where the equivalent width seems to be independent of the luminosity. Although the data has 
a large dispersion, part of it is probably due to intrinsic variability. Plotting the time average 
equivalent width vs continuum luminosity gives a significantly smaller dispersion (Kinney, Rivolo 
and Koratkar 1990, fig 4). From these data we get EWOciv = 120 A. 



The slope of the Baldwin effect predicted by eq. |37| is 



"'"'f!"'-' = -lOE-L^ (38) 

(llnLt 46 

For E c = 20eV this gives a slope of ~ O.SL^g 4 . 

Fig. 2 shows the data referenced above, along with the model prediction (eq. |37j) as a dashed 
curve. We note that the decrease in the equivalent width predicted by eq. [37] is too large for 
luminous quasars. 
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The Baldwin Effect - Data and Models 
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Figure 2. The Baldwin effect (equivalent width of CIV vs. the monochromatic continuum 
luminosity at 1550A. The fuzzy thick curves indicate the dispersion of the data, taken from 
Kinney, Rivolo and Koratkar (1990). The smooth curves indicate the prediction of the models 
derived in this work: heavy dashed - a power law ionizing spectrum, L u ~ dotted - with 

constant L/M (eq. |37] with E c = 20 eV) , solid with L/M oc L 1 / 2 (eq. ||). 

This is a result of our assumption that L/M ~ constant. Actually, as discussed in section 4 
the L/M ratio appears to be correlated with luminosity (see fig. 1), with the exception of NLS1. 
A similar result is obtained by Wandel and Petrosian (1988) and (although on a small range in 
luminosity) by Wandel, Peterson and Malkan (1999, fig. 4). 
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From fig. 1 we adopt an approximate relation 



L/L Edd « 0.1L 



1/2 



16 



Combining this with eqs. yields 



1/8 



E m (20 eV)L 



16 



which gives for the equivalent width 



EWciv = EW0 C ivexp{-2L l 4 <°). 



(39) 



This result is shown by the solid curve in fig. 2, which is consistent with the observed Baldwin 
effect . 

For Lya we have Ej — Ei on = 9.2 — 13.6 = 4.4eV, which gives 

* w (f*V = -I.IE^L%\ (4( 

which is consistent with the weaker Baldwin effect found for the Lya line (Kinney, Rivolo and 
Koratkar 1990). 



8.1. Comparison with the Data: Luminosity and Energy dependence 

Our derivation predicts a weak luminosity dependence of the Baldwin effect . While this 
dependence (eq. |3^) on luminosity, indicating a slow increase of the Baldwin effect (or a decrease 
of the slope of the line luminosity vs. continuum luminosity) for very luminous quasars is not 
obvious from the data, it cannot be ruled out either, because of the large scatter. When the 
scatter is reduced by averaging multiple observations of the same object, (e.g. fig. 4a in Kinney, 
Rivolo and Koratkar (1990) the resultant Baldwin effect seems to suggest an increase the slope of 
EW plotted vs. L at high luminosities. 

We have seen however, that the increase in the Baldwin effect predicted by eq. |38| is too 
strong at the highest luminosities. When the increase of the Eddington ratio with luminosity is 
taken into account, the model can fit the CIV data for the highest energies, for E c ~ 20eV . This 
value may be obtained for a disk extending to the last stable orbit or by including Comptonization. 

Alternatively, the excess effect could be compensated by the fact that the parameters m 
and tl may depend on luminosity; If fn increases with L (as indicated in Wandel, Peterson and 
Malkan 1999) or the variability timescale increase with L less than linearly, E c will increase with 
luminosity, compensating for the weak luminosity increase. An additional uncertainty is due to 
the unknown dependence of the emissivity (which we assume to be constant) on L, which may be 
influenced, among other factors, by the covering factor. 



8. Discussion 
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8.2. Comparison with the Data: the Baldwin effect in other lines 

The Baldwin effect has been observed also for other lines. Our model predict a strong 
dependence of the Baldwin effect on the ionization energy. It is interesting to see whether this 
prediction could reproduce observed strength of the Baldwin effect in different lines. Fig. 4b of 
Kinney, Rivolo and Koratkar (1990) shows a much weaker effect for Lya, and in particular for 
low luminosity objects (L > 10 4 6erg/sec) the Lya data are consistent with a constant equivalent 



width , which is predicted by our model (eq. 4C). 



By the same token, our model predicts that lines with a higher ionization potential like OVI) 
should have a stronger Baldwin effect , which is also observed (Korista et aZ.1995). 

As the emissivity of CIV and OVI has a similar characteristic distribution in the Q — n plane, 
one could expect a similar Baldwin effect for the two lines, which is indeed observed. The lack 
of that sharp feature in the emissivity of lines such as Lya and H/3 (Korista et al. 1998a) could 
explain the fact that these lines do not show a significant Baldwin effect or only a weak one. 
Such a specific functional form can be predicted by assuming an explicit radial distribution of the 
density and covering factor. This the expected functional form of the Baldwin effect could provide 
an observational test of the main suggestion of the paper. 



8.3. Normalization 



We have reproduced the Baldwin effect from the basic relations between black hole mass, 
accretion disk temperature combined with the line emissivity. It would be of interest to estimate 
the normalization for the proportionality relationship found, which would enable a quantitative 
prediction of the expected Baldwin effect . However, the relation between the emission line 
luminosity and the ionizing continuum shape depend on the unknown density profile of the line 
emitting gas, and on the unknown differential covering factor, as discussed in section 5. In 
the absence of more detailed data on the radial distribution of the line emitting gas, further 
analyses is limited to qualitative predictions, for which the positive correlation with the radiation 
temperature demonstrated in sections 5-7 is sufficient. A more explicit dependence predicted by 
photoionization models can be obtained using the equations derived in sections 5 and 6, if the 
radial distribution (characterized by n(r) and df(r)/dr) are known. 



9. Summary 

We suggest that the Baldwin effect is a result of the distribution of black hole masses 
combined with the thermal nature of the line-driving ionizing continuum. Using the variability- 
black body relation, or alternatively the thin accretion disk spectrum, we derive an analytical 
estimate for the dependence of the continuum spectrum on the black hole mass and show that the 
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ratio of ionizing continuum to the continuum luminosity near the lines (related to the equivalent 
width ) is slowly decreasing with increasing mass. We find that the for a model- independent 
black body temperature and for most variants of the thin accretion disk model the equivalent 
width decreases with increasing continuum luminosity, reproducing the observed Baldwin effect . 

Stimulating discussions with Jack Baldwin during the meeting "Quasars as Standard Candles 
in Cosmology" at La Serena, Chile, triggered this work. I also acknowledge valuable discussions 
with Gary Ferland, Matt Malkan and Vahe Petrosian, the contribution of an anonymous referee, 
and the hospitality of the Astronomy Department at UCLA and Stanford University. 
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